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ABSTRACT 


Field measurements from a cross~—shore array of nine pressure 
sensors, spanning the surf zone, are used to examine the evolution of 
ensemble averaged wave face slopes of ocean waves as they propagate 
through the breaking region. Averaged wave slopes are determined from 
time series of the measured sea surface elevation and from an averaged 
waveform calculated from bispectral coefficients, and compared with 
predictions from a wave transformation model that includes wave breaking 
described by rollers. Measured percent wave breaking are used to 
examine the evolution of third moment and bispectral statistics in 
relation to breaking patterns. 

Shoaling waves gradually transform from peaked, Stokes-like waves 
to forward pitched asymmetric waves just prior to breaking. Inside the 
Surf zone, wave asymmetry is modified by the breaking distributions and 
the effects of bottom topography. The observations Suggest a relation- 
ship between the cross-shore wave breaking distributions and wave 
Slopes. Wave slopes predicted using a calibrated wave transformation 
model which includes wave rollers are in qualitative agreement with 


measured wave slopes. 
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I. INTRODUCTION 


AS waves shoal in intermediate water depths, nonlinear 
interactions transfer incident wave energy from the peak 
incident frequency to both higher (harmonics) and lower 
(infragravity) frequencies. The strength of these interac- 
tions depend on many factors, which include incident wave 
properties, topography and whether or not breaking is 
occurring. On the inner shelf, wave energy flux is approxi- 
mately conserved. However, in very shallow water, shoaling 
waves become unstable and eventually break (creating the 
Surf zone), and organized wave energy is converted to 
turbulence and dissipated. Turbulence generated in the 
Surface boundary layer by wave breaking is the primary 
dissipative mechanism in the surf zone. To model wave 
transformation all the way to the shoreline requires that 
the dissipative mechanism be specified. 

A common approach to describing breaking wave dissipa- 
tion has been to use the analogy of a hydraulic jump to 
bores in the surf zone (LeMehaute, 1962; Engelund, 1981). 
This approach has been successful at predicting wave 
transformation across arbitrary bathymetry (e.g., Battjes 
and Janssen, 1978; Thornton and Guza, 1983; Roelvink, 1993; 
Lippmann, et al., 1996). Many of these models use empirical 
weighting functions to describe the portion of wave break- 


ing. Although they do a good job of predicting wave 








heights, the spatial distribution of wave breaking is not 
always well represented. an improvement to the wave 
transformation model using the concept of wave rollers (as 
described in Duncan, 1981; Engelund, 1981; Svendsen, 1984; 
and others) has resulted in a better representation of 
breaking distributions, as shown in Figure 1 by comparing 
the roller model by Lippmann and Thornton (1996), henceforth 
referred to as LT96, (solid line) with a typical wave 
transformation model (dashed line). The formulation of wave 
breaking in roller models (Svendsen, 1984; Nairn, et al., 
1990; Dally and Brown, 1995; LT96; and others) suggests that 
the slope of the front face of the wave is important in 
determining the spatial variation in wave breaking distribu- 
tions in the surf zone. 

These models incorporate roller energy advection terms 
in an energy balance equation that retains the energy flux 
within a turbulent mass of water perched on the wave face 
(the roller) as waves break, advects roller energy landward, 
and then loses energy to turbulent dissipation in the 
underlying water column by action of the shear stresses at 
the wave/roller interface. The energy flux balance is 


described by 


wid + eee ae = “TSC , (1) 


where EF, and E. are wave and roller energy, C, is the group 


2 




















velocity, C is phase speed of shallow water waves, and t is 
the shear stress at the wave/roller interface. Both the 
roller energy and shear stress depend on the roller area, A 
(Duncan, 1981; Svendsen, 1984; Deigaard, 1993; and others). 
The size of A is determined from the energy lost in a 
hydraulic jump of height, H. The area is determined by the 
energy lost to turbulence while conserving mass and momentum 
both upstream and downstream of the flow discontinuity 
(Engelund, 1981). The area per unit crest length of the 
roller describing the turbulence generated in the jump is 


given by 


3 
eccre 7 
where h is the water depth and o is the angle of jump 
discontinuity. The spatial variation of the rollers thus 
depends on o. A schematic showing the geometry of a wave 
with a turbulent roller on a linear wave face is depicted in 
Figure 2. Using eppie geometrical approximations (Figure 


2), 0 18 given by a function of breaking wave steepness 


(HL/L) 


ene) ey (3) 
L 


where ( ) indicates ensemble averaging, L is the local 


wavelength, and ” (which is equivalent to B/e where B and e 





are the fractions of wave height and wave length covered by 
the roller) is an assumed spatially invariant constant 
dependent on how skewed or pitched forward the waveform. 

The value of yw is a free parameter and used as a 
calibrating constant to give a best fit to the observations 
of concurrent sea surface elevation and wave breaking 
statistics (obtained by pressure sensors and video record- 
ings of the surf zone). The wave data is used to estimate 
how well Equation 3 resembles the actual breaking wave 
slopes. 

The objective of this paper is to examine evolution of 
the wave face slope across the breaking region and compare 
the observations with the LT96 roller model predictions to 
assess the formulation of the model. Two methods, time 
series and bispectral, are used to calculate wave slopes 
from pressure sensor data. Slopes calculated directly from 
the time series include contributions by all incident waves. 
In the bispectral method, a slope is calculated from a 
Single averaged waveform composed of a primary wave and its 
nonlinearly coupled harmonics whose amplitudes and phases 
are determined from the bispectrum. The roller model 
predicted wave face slopes are compared to both the average 
Slopes calculated from time series and bispectral methods. 

The next section describes how statistics from pressure 
sensors are used to calculate wave face slopes. In the 


experiment section, the field measurements of sea surface 














elevation and wave breaking obtained during the 1990 Delilah 
experiment are described. In the results and discussion, 
the evolution of the qualitative behavior of wave 
nonlinearities are examined as waves propagate shoreward 
through the surf zone. Observations of percent wave 
breaking, nonlinear third moments (skewness and asymmetry), 
and bispectral statistics are examined in relationship to 
bathymetry. The cross-shore behavior of ensemble averaged 
wave slopes is compared with predictions from the LT96 


roller model calibrated with observations of wave breaking. 











II. ANALYSIS 


AS waves shoal in intermediate water depth, they become 
Stokes-like with peaked crests and lengthened, shallow 
troughs. As waves enter progressively shallower depths, 
they begin to pitch forward until they break creating the 
Surf zone. Across the Shoaling and breaking regions, energy 
at the primary frequencies is transferred nonlinearly to 
both lower and higher frequencies. The higher frequency 
forced waves are phase-coupled to the primary. The degree 
of coupling and phase relationship between primary and 
harmonic frequencies result in wave deformation leading to 
breaking and are quantified with higher (third) order 
Spectral techniques (the bispectrum) . 

The waveform evolution of Shoaling surface gravity 
waves is described using the third moment wave Statistics, 
skewness and asymmetry, of the time series (Elgar and Guza, 
1985). Profiles for both skewed and asymmetric waveforms 
are shown in Figure 3 (taken with slight modification from 
Doering, 1988). Skewed waves are similar to Stokes waves, 
characterized by a peaked wave profile with elongated and 
flattened troughs, and a lack of symmetry about the horizon- 
tal axis. Asymmetric waves resemble a pitched forward 
profile like "saw-tooth" bores common to the surf zone. 

Both skewed and asymmetric waves have the same amplitude 


Spectrum as illustrated by the Fourier decomposition in 








Figure 3. The phase shifting of harmonic components 
relative to the primary frequency is a characteristic of 
asymmetric and skewed waves. The difference is in the phase 
Spectrum, also shown in Figure 3. The skewed waves have 
phase-locked harmonics in-phase with the primary frequency. 
A purely asymmetric waveform, or saw-toothed profile, has 
the n* harmonic phase shifted 1/2 relative to the primary 
frequency. All other phase relationships result in waves 
having both skewed and asymmetric characteristics. 

First order spectral techniques cannot be used to fully 
examine these higher order moments. Higher order spectra 
(bispectrum) have been applied to ocean surface gravity 
waves to study nonlinear phenomena for several decades 
(Hasselman, et al., 1963, and others). Elgar and Guza 
(1985) calculated bispectra of shoaling waves and found that 
harmonic growth and biphase evolution across the shoaling 
region was consistent with gross changes in wave shape, and 
used this technique to study the evolution of both skewness 
and asymmetry in the shoaling region. 

If we represent waves entering the surf zone by a 
stationary Gaussian function of time with zero mean, &(t), 
then the wave field is represented by a superposition of 
Statistically independent wave components. Thus, the time 
series is completely described by the double-sided energy- 


density spectrum, S(f), expressed in terms of Fourier 


coefficients, F(f£), 

















S(f) = EL[F(f) F'(f)] / af (4) 


where f is frequency, Af is the frequency bandwidth, and 
E{] 1s the expected value operator. The variance of the 
time series is simply the integral of the energy-density 


spectrum 


co 


E[&(t)] = [s(£)df (5) 


-® 


The bispectrum can also be expressed in terms of 
Fourier coefficients (Hasselmann, et al., 1963; Kim and 


Powers, 1979) 


B(Y,7f) = 29) FS) S (£485). 7 A (6) 


The bispectrum is nonzero only if waves at frequencies f, 
and f, are nonlinearly coupled, generating both sum and 


difference frequencies, f,, 
£2 2,24, (7) 


Unlike the energy-density spectrum which is real only, 
the bispectrum is complex, having real and imaginary 


components which are related to third moment statistics. 











The integral of the imaginary part (S[]) of the bispectrum 
normalized by the time series variance to the 3/2 power, is 


used to define the asymmetry 


Oo oo 


[ [StB(4,,f,) ld£,aE, (8) 


Asymmetry = ELE (te) 


equivalent to the skewness of the Hilbert Transform (see 
Bendat and Piersol, 1986) of the time series (Elgar and 
Guza, 1985). Similarly, the normalized integral of the real 


part (3[]) of the bispectrum is used to define the skewness 


co 2 


[ [ R1B(4,,£,) 1d£,aF, (9) 


Bete) ee 


Skewness = 


equivalent to the third moment of the time series. 
Bicoherence, b, is the normalized form of the 
bispectrum (0 < b < 1) and is used to Separate independent, 


free waves from the nonlinearly phase-coupled waves (Kim and 


Powers, 1979) 


b?(£,,f,) = mad 13 


(10) 
E[\7(£,) F(£,)/°] BUF (f,+£,))7) 
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Values of bicoherence greater than the significance level 


for a given confidence interval indicates the presence of 
nonlinear coupling between interacting frequencies. The 
Significance level for 95% confidence on nonzero 


bicoherence, b,, is calculated as 


5? 


6 


Be. S 
°° ‘dof 


(11) 
where dof is the number of degrees of freedom (Elgar and 
Guza, 1985). 

The biphase represents the relative phase relationship 


between the interacting frequencies and is determined by 


ce 
B(f,,f2) te ee ei} (12) 


H[B(f,,f,)] 


For significant bicoherent coupling at (f,, f,), biphase 
values of 0 represent contributions to a more skewed 
waveform (Figure 3), whereas values of -90° make the 
waveform more asymmetric (Figure 3). Other values of 
biphase indicate harmonics are contributing to both skewed 
and asymmetric waveforms. 

Surface elevation slopes, 07/dx, are estimated from the 


time series of the sea surface elevation, y(t), by 


li 





on _ 10n ‘3 
c ot oa 


where x 1S cross-shore coordinate and t is time, and 


assuming that the wave field is non-dispersive with phase 
speeds dependent only on the (known) local depth, c = Vgh 


with g gravity. The differentiated surface elevation time 
series, dn/dt, is obtained in the frequency domain by 
applying the transfer function of -iw to the complex 
Fourier amplitudes of y(t) and inverse transforming back to 
the time domain. The surface slope is calculated as the 
maximum wave face slope found from the extrema of a°n /at?. 


The second derivative of the surface elevation time series 
is obtained by applying the transfer function, i2w?, to the 


complex Fourier amplitudes of y(t) and inverse transforming. 
The up-crossings are identified and the maximum positive 
Slopes between the up-crossings are the slopes of the wave 
face. The ensemble averaged wave face slope is calculated 
as the rms slope. 

AS a second technique, the wave face slope is estimated 
from the averaged waveform represented as a primary wave at 
the peak incident frequency plus the nonlinearly coupled 


harmonics found from the bispectrun, 





(€) = €2.coS(w .t) + a,b,cos(2e,t - $,) + a,b,cos(30,t ~ $,) + ..(14) 
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The time series represented by Equation 14 is periodic with 


frequency f£,. Phase locking makes all the waves appear to 
be the same. The same procedure for calculating the wave 
Slopes used in the time series analysis is applied to 
Equation 14 over a single wave period. Thus, the maximum 
0n/OX Over one wave period defines the average angle of the 
wave face. The wave slope is now a function of a primary 
(aj) and its coupled harmonic (Q,, a, -..) amplitudes, 
coupling coefficients or bicoherences (b,, bj, ...), and 


biphases (4, Poy---)- 
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III. FIELD EXPERIMENT 


The data were collected as part of the Delilah 
Nearshore Processes experiment held at the U. S. Army Corps 
of Engineers Coastal Engineering Research Center Field 
Research Facility (FRF) in Duck, North Carolina during 
October 1990. Nine pressure sensors were deployed in a 
cross-~shore array from ~4.5m depth to the shoreline (Figure 
4). Data were acquired approximately continuously at 8 Hz 
for about 3 weeks. 

Bottom pressure measurements were converted to sea 
surface elevation by Fourier transforming the one hour 
pressure record, applying the linear wave theory spectral 
transformation function, and inverse transforming to obtain 
sea surface elevation (Guza and Thornton, 1980). At the 
same time, the data were band-pass filtered from 0.05 to 0.5 
Hz by zeroing the Fourier coefficients prior to inverse 
transforming to remove high frequency noise and lower 
frequency (infragravity) waves. The sea surface elevations 
were detrended, and subsampled to 2 Hz for subsequent 
analysis. 

Video recordings of the surf zone were obtained from 
cameras mounted on a 44m high observation tower. The cross- 
Shore distribution of the number of waves breaking at a 
particular location were estimated using the video taken 


during daylight hours. The number of waves breaking were 


15 





counted manually from video pixel time series (Lippmann and 
Holman, 1991) at approximately 10m intervals Spanning the 
width of the surf zone. 

The bottom profile was surveyed daily using the Coastal 
Research Amphibious Buggy (CRAB; Birkemeier and Mason, 
1984). The bathymetry evolved Significantly during the 
experiment in response to the variable wave climate, which 
included two events exceeding 2m offshore Significant wave 
height, H,. A single, well defined sand bar migrated 
Seaward as the waves increased. Offshore wave conditions 
were generally narrowband in both frequency and direction 
(measured in 8m depth with a linear array of bottom mounted 
pressure sensors). For a more complete description of the | 
Delilah experiment, see Birkemeier, et al. (1991) and | 
Thornton and Kim (1993). 

Two cases are examined in detail to contrast wave 
transformation over a monotonically increasing profile with 
a barred profile. Data to be presented are from a period 
when bathymetry was in a terrace-like configuration (8 
October 1212 Hrs. at high tide) during which the wave field 


was narrowband (f, = 0.093 Hz) with H, * 0.51m; and during a 





period when the bottom was in well-developed bar-trough 
sequence (12 October 0718 Hrs. at low tide) with a more 
energetic (H, ~ 0.93m) and more broadband (0.09 < £f, < 0.12) 
wave field. All spectra and bispectra are computed with 112 


degrees of freedom. 
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IV. OBSERVATIONS 


The sea surface elevation spectrum for 08 October, 
taken at the most offshore sensor located in about 4.5m 
depth (Figure 4) shows the waves were relatively narrowband 
(Figure 5b). Contours of bicoherence greater than the 953 
Significant level (Bs: = 0.054 calculated as in Equation 11 
for 112 d.o.f.) plotted in bi-frequency space are shown in 
Figure 5a. Only the unique portion of the bispectral domain 
is shown (Kim and Powers, 1979). The high bicoherence at 
the peak of the spectrum (f. 0.093 Hz) indicates there is 
Significant nonlinear generation of the first harmonic (at 
frequency 2f,) by self-self interaction of the primary. The 
energy spectrum for 12 October (Figure 6b) is broader than 
the spectrum for the 8th. The nonlinear interactions are 
Spread out over a range of incident frequencies, as shown by 
the bicoherence contours in Figure 6a. 

The cross-shore evolution of the normalized third 
moment statistics (wave skewness and negative asymmetry) 
across the bottom profile are shown in the lower panel of 
Figure 7 for 08 October. The negative asymmetry is due to 
the selection of the positive offshore coordinate reference 
frame. The observed spatial distribution of the percent 
wave breaking is shown in the upper panel. The observed 
breaking region starts from about 80m offshore and increases 


monotonically towards the shoreline where nearly all the 
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waves are breaking. The skewness is high offshore and 
increases Slightly as the narrowband waves propagate 
Shoreward. On the other hand, the asymmetry is nearly zero 
offshore, but increases rapidly beyond the break point as 
waves shoal, suggesting the waves are pitching forward just 
prior to breaking (similar to results obtained by Elgar and 
Guza, 1985). Inside the break point, the asymmetry remains 
high all the way to the shoreline. 

Percent breaking, skewness and asymmetry for the well- 
developed barred profile (12 October) are shown in Figure 8. 
Waves began to break about 130 meters offshore, increasing 
to a maximum of over 60 percent at the crest of the bar. 
The percentage of wave breaking decreases as waves propagate 
into the deeper water of trough of the bar, reaching a 
minimum of less than 5% at the Shoreward limits of the 
trough and then increases again as the waves break on the 
foreshore. 

The skewness is high offshore, increasing gradually as 
waves shoal, break on the bar, and reform in the trough. 
Only very near the shore-break does the skewness decrease. 
The asymmetry is small offshore (Similar to the 8th) and 
again increases Sharply just prior to breaking. The 
asymmetry is largest over the bar crest where the maximum 
percent of waves are breaking. Inside the bar crest, the 
asymmetry increases again at the shore-break. 


The cross-shore changes in the relative amplitudes of 
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the primary, a,, occurring at frequency foe and the first 
two harmonics, a,b, and a,b,, (obtained from Figures 5a and 
6a), occurring at frequencies 2f,, and 3f,, are shown in 
Figure 9 for 08 October (plotted as the ratios a,b,/a, and 
a,b,/a,)- As the waves propagate shoreward the amplitudes of 
the harmonics grow relative to the primary, increasing more 
rapidly inside the surf zone where breaking limits the size 
of a,;. The bicoherences and biphases between the primary 
and harmonics are also shown in Figure 9. The bicoherences 
are strong from offshore through the surf zone showing the 
wave field to be significantly coupled nonlinearly. The 
relative phase of the harmonics to the primary (the biphase) 
is nearly zero offshore (in-phase) where the waves are 
peaked up but not pitched forward. The biphase of both 
harmonics increases gradually as waves approach the surf 
zone, continues to steepen as waves begin to break, and 
approaches a maximum biphase of about -7/2 in the inner surf 
zone. 

The amplitude ratios, bicoherences, and biphases for 
the broader, more energetic (12 October) data are shown in 
Figure 10. The harmonic amplitudes are high both offshore 
and in the surf zone, including the trough region shoreward 
of the well-developed bar. The bicoherences are also high 
throughout the shoaling and breaking regions, except very 
near the shoreline. The biphases of the harmonics again 


increase from approximately zero offshore to near -1/2 over 


ig 








the bar crest. The change toward smaller biphase of the 2nd 
harmonic, while its amplitude and coherence remain high, 
Suggests significant nonlinear evolution in the trough of 
the bar, leading to a less pitched forward wave. 

The root mean square (rms) wave face Slopes calculated 
from the derivatives of the time series and from the maximum 
Slope determined from the averaged waveform calculated using 
the bispectral coefficient (Figures 9 and 10) are shown in 
Figures 11 and 12. For 08 October (Figure 11), both methods 
Show an increase in wave face angle as waves approach and 
then enter the surf zone. Offshore the estimates from the 
time series are lower than from the bispectral method, not 
Surprising considering that time series includes all waves 
whereas the waveform expressed by Equation 14 includes only 
phase-locked harmonics and excludes other coupled wave 
motions. Wave slopes are maximum near the outer limits of 
the surf zone, and decrease slightly as wave amplitudes are 
reduced by breaking. The wave slope behavior for 12 October 
is similar to 08 October, with maximum wave Slope near the 
bar at the outer edge of the surf zone (Figure 12). 

The predicted wave slopes computed from the roller 
model (Equation 3) calibrated to best fit observations of 
percent breaking are also shown in Figures 11 and 12. The 
upper panels show the best fit for yw = 4.35 and wy = 4.05 
for 08 and 12 October. The value of H, in Equation 3 is 


obtained by multiplying the H_. = v8 o,, where o, is the 
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surface elevation variance, by the Rayleigh distribution and 
a spatial weighting function calibrated from percent 
breaking observations (Lippmann and Thornton, 1996, Equation 
23). The wavelength (L) is computed from £, and local depth 
using the shallow water approximation. In the surf zone the 
cross-shore behavior of the roller model predicted slopes is 
in qualitative agreement with the observed slopes. The 
roller model slopes are maximum in the outer limits of the 
surf zone, and decrease in the trough. The qualitative 
agreement suggests the roller model does a reasonable job of 


representing wave breaking in the surf zone. 
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V. DISCUSSION 


Wave slope predictions are dependent on the value of 
wy, which in turn is dependent on the accuracy of the wave 
and roller transformation model. This means that accurate 
estimates of wave energies and percent breaking distribu- 
tions are vital steps in accurately predicting wave slopes. 
The presence of the bar tends to complicate the predictions 
by causing a region of intense breaking followed by a region 
of decreased breaking. The spatial misrepresentations in 
the breaking region alter the yw value which gives the best 
fit to the percent breaking observations (Figures lla and 
12a). Fortunately wave heights and breaking observations 
are usually well modeled using the energy balance equations. 
Thus, differences between calculations from observations and 
predictions in the estimated value of yw are minor (Lippmann 
and Thornton, 1996). The effects on wave slope predictions 
are also minor. 

The LT96 roller model assumes constant yw across the 
bottom profile when estimating the wave slopes. In their 
model, wave face slopes are a function of the local, 
ensemble averaged breaking wave steepness (Equation 3). The 
data and model estimated wave face slopes are in general 
qualitative agreement (particularly in the surf zone), 
suggesting that the wave face angles are indeed proportional 


to wave steepness. The changes in steepness are the result 
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of the cumulative effects of both the evolution of the wave 
profile from a skewed, Stokes-like waveform to a more 
pitched forward, saw-toothed like profile, and from the 
dissipative effects of wave breaking which severely limit 
wave heights, but not wavelengths, in the surf zone. 

The roller model predicts the wave face Slope of 
breaking waves only, whereas the time series measurements 
include all waves, pusdiine and non-breaking. The 
bispectral waveform describes only the shape of the primary 
waveform, and thus should be better representative of the 
breaking waves wave slopes. Therefore, outside the surf 
zone, model and bispectral measured wave Slopes are expected 
to be greater than the averaged wave face slopes from all 
waves of the time series (Figures 11b and 12b). Inside the 
Surf zone where most of the waves are breaking, the roller 
model predictions are in good qualitative agreement with the 
root mean Square (rms) averaged wave face Slopes. The wave 
face slope comparisons are reasonable for both narrowband 
waves over a monotonically increasing bathymetry and the 
more complicated shoaling and breaking over a barred beach. 

Both measured and roller modeled wave face Slopes 
(Figures 11b and 12b), which range about 4-10 degrees (6 
degrees is about a 1:10 Slope for comparison), seem small 
compared with casual observations of breaking waves which 
appear to approach vertical for plunging breakers. The very 


Steep wave faces observed on large breakers are horizontally 
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small in dimension compared with the wavelength. Thus 


sampling the very narrow, highly nonlinear wave crest is 
important, and has proved limiting using standard wave 
measuring techniques (i.e., bottom pressure meters), which 
depend on the sensor resolution and sampling rates (typi- 
cally 2-8 Hz for surf zone applications). Even at much 
higher rates it is still not clear that local time series 
are adequate, owing to hydrodynamic filtering of the very 
high frequencies (which are attenuated even in shallow 
water) that produce the steepness in the upper crest of the 
wave. 

Direct measurements of wave slopes in the crest of the 
wave requires spatially resolving measurement techniques, 
which owing to logistical difficulties in deploying a large 
array of fixed instruments, will most likely have to have a 
remote sensing nature (such as stereo video). The direct 
Spatial measurements also eliminate the need to estimate the 
local phase speed necessary when using single point time 
series (Equation 13). 

Ensemble averaged phase speeds in shallow water are 
modeled with linear theory (Thornton and Guza, 1982; 
Lippmann and Holman, 1991; and others), which tend to 
underestimate actual phase speeds due to some amplitude 
dispersion. Underestimation of phase speed results in 
overestimating the wave face slope magnitude using Equation 


13. Including amplitude dispersion would greatly complicate 
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the analysis. The linear approximation, however, is 
sufficient for observing qualitative wave face slope 
changes. 

Wave face slopes change as a result of evolution of 
skewness, asymmetry, and wave steepness, which in turn are 
determined by the shoaling topography (Figures lic and 12c) 
and wave breaking distributions, which may be highly 
variable in the surf zone. Observations of wave asymmetries 
and percent breaking (Figures 7b and 8b) appear to have 
Similar spatial distributions (Figures 7a and 8a). Even 
though wave breaking is expected to Significantly modify the 
wave steepness, the exact effect on wave skewness and 
asymmetries can not be resolved. 

Bispectral analysis leads to estimates of wave slopes 
through estimates of harmonic amplitudes (Figures 9b and 
10b), and the degree of coupling between waves and their 
phase relationships (Figures 9c and 10c) at interacting 
frequencies. Estimates of these quantities can be sensitive 
to the method of computing the bispectrum (Elgar and Guza, 
1988). In this work, the filtered time series are broken up 
into 512 point samples (corresponding to 256 second re- 
cords), a Hanning data window applied, and then are averaged 
with no frequency merging done. The effects of frequency 
bandwidth on wave slope calculations derived from bispectral 
estimates was tested by using a range of ensemble widths 


from 256 to 2048 points (corresponding to 448 to 56 degrees 
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of freedom). The resulting deviation in wave slope was 


slight, with qualitative trends across-shore being pre- 
served, and about a +0.5-1.0 degree variation occurring 
occasionally about an average wave slope which was not 
substantially effected. Varying the high (and low) fre- 
quency cutoff’s from the band-pass filtered applied to the 
raw time series did not affect our estimate of slope in any 


of our methods. 
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VI. CONCLUSIONS 


Ensemble averaged wave face slopes of breaking waves in 
the surf zone are calculated directly from the time series 
and from an averaged waveform calculated from bispectral 
coefficients and compared with predictions from a wave 
transformation model that includes wave breaking described 
by rollers. Also examined is the relationship between 
percent breaking and both normalized third moments (skewness 
and asymmetry) and bispectral statistics. 

The nonlinear evolution of waveform across the breaking 
region is qualitatively described from observations of third 
moment statistics and from examination of the phase rela- 
tionships between the primary frequency and its coupled 
harmonic components. The roller model predictions are in 
qualitative agreement with the observed waveform calculated 
from observations. 

The skewness is generally high across the shoaling and 
breaking regions, consistent with a peaked waveform common 
to Stoke-like gravity waves. As the waves approach breaking, 
the wave crests become more pitched forward, transforming 
from near zero asymmetry offshore (indicating little or no 
pitch) to a maximum in the inner surf zone (indicating a 
more pitched forward waveform). The asymmetry reaches a 
maximum coincident with the location of breaking. 


Bispectral analysis shows that coupled harmonic amplitudes 
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grow through the shoaling and breaking regions relative to 
the amplitude of the primary incident frequency. Shoreward 
of a well-developed bar, the waves gradually reform from 
breaking pitched forward bores to skewed non-breaking waves. 
It is not clear that this transition is due to topographical 
influence or wave breaking distributions, or both to some 
degree. 

For the barred profile the average wave face angle 
first increases as waves pitch forward over the bar just 
prior to breaking, then decreases as waves propagate into 
the trough, and then increases again as waves break at the 
shoreline. The wave slopes predicted from a calibrated 
(with percent breaking observations) roller transformation 
model are in qualitative agreement with the wave slope 
observations, suggesting that roller theory well represents 


waveform evolution particularly in the breaking region. 


30 





APPENDIX 


Figure 1. Roller model predictions of percent wave breaking 
from video recordings of the surf zone from 10 October, 1990 
during the Delilah experiment (From Lippmann and Thornton, 
1996). The upper panel shows percent breakers observed in 
the field (circles) and predicted by the roller model (solid 
line) on the nearshore profile shown in the lower panel. 
Also shown for comparison are predictions from a (typical) 
previous wave transformation model (dashed line, Lippmann, 


et al., 1996). 


Figure 2. Schematic showing geometry for a wave with 
turbulent roller on the linear wave face (modified slightly 
from Lippmann and Thornton, 1996). Variables are defined in 


the text. 


Figure 3. Schematic (modified from Doering, 1988) showing 
Sinusoidal (horizontal and vertical symmetry), skewed 
(horizontal symmetry; vertical asymmetry) and asymmetric 
(horizontal asymmetry; vertical symmetry) waveforms. 
Amplitude and phase relationships relative to the primary 
frequency are also shown to the right of each wave profile. 
Note that although the amplitude spectra are identical for 
the skewed and asymmetric waves, the harmonics for the 
skewed wave are in-phase, whereas the harmonics are 1/2 


phase shifted relative to the primary for the asymmetric 
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wave. 


Phase relationships between primary and harmonics 


which are between 0 and m/2 have both skewed and asymmetric 


characteristics. 


Figure 4. Beach profiles for 08 and 12 


October during the 


Delilah experiment and the location of pressure sensors. 


Cross-shore distance is relative to the 


Field Research 


Facility (FRF) coordinate System and depths are relative to 


National Geodetic Vertical Datum (NGVD). 


Figure 5. 

for 08 October at 1212 Hrs. Observed at 
sensor (Figure 4). 

Series with 56 d.o.f. 


Of 0.0075 Hz. 
(greater than the 952 Significant level 
bi-frequency space corresponding to the 


Spectrum in (a). Contour interval is 0. 


bicoherence of 0.79 occurring at f, £ 


2 


Figure 6. (a) Same as Figure 5a for 12 


(b) Same as 5b with maximum bicoherence 


Seal 


, = 0.094 Hz and f. 


£ 0.125 Hz. 


Figure 7. Cross-shore evolution of the 
moments skewness and negative asymmetry 


08 October. Statistics are computed in 
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(a) Sea surface elevation energy-density spectrum 


the most seaward 


Spectrum is computed from one hour time 
corresponding to a frequency bandwidth 


(b) Contours of significant bicoherence 


of 0.23) plotted in 
energy-density 


1 with maximum 


0.094 Hz. 


October at 0718 hrs. 


of 0.48 occurring at 


normalized third 
(center panel) for 


the time domain from 


ee 




















bandpass filtered sea surface elevation time series. The 
bathymetry is shown in the lower panel and the observed 


percent wave breaking in the upper panel for comparison. 


Figure 8. Same as Figure 7 for 12 October at 0718 hrs. 


Figure 9. Cross-shore evolution of (b) observed ratio of 
primary to first (Stars) and second (circles) harmonic 
amplitudes, (c) negative biphase between primary and first 
(stars) and second (circles) harmonics, and (d) bicoherence 
between interacting primary frequencies (stars) and interac- 
tions between primary and harmonic frequencies (circles) for 
08 October. The bathymetry is shown in (e), and the 


observed percent breaking in (a) for comparison. 


Figure 10. Same as Figure 9 for 12 October. 


Figure 11. Cross-shore evolution in rms wave face slopes 
(middle panel) calculated from the time series and maximum 
Slope from nonlinearly coupled Equation 14 for 08 October. 
Observed (stars) and roller model predicted (solid line) 
percent breaking is shown in the upper panel. The roller 
model is from Lippmann and Thornton (1996) and has best fit 
y parameter if 4.35. Also shown in the middle panel is the 
roller model generated wave slopes based on the yw value 


obtained (Equation 3). The bathymetry is shown in the lower 
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panel for reference. 


Figure 12. Same as Figure 11 for 12 October. Best fit y 


value from roller modeling is 4.05. 
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